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Surface  atomic  structures  of  clean,  oxygen-adsorbed .  and  oxidized  (111)  nickel 
have  been  stuuiej  quantitatively  bv  using  MeV  ion  scattering  in  eombinati'a 
with  AF.S  and  NEED.  We  show  that;  the  clean  (111)  nickel  surface  has  the  bulk- 
like  st rue  fare  with  re  eons t rue t ion  or  relaxation  less  than  0.02  A,  the  surface 
thermal  vibration  amplitude  is  enhanced  by  'MO;;  compared  to  the  bulk  value, 
adsorbed  oxvveu  results  in  surface  lattice  expansion  by  Ml. 15  A  wiiicu  is  closet 
correlated  tv~  the  p(2::2)  or  ( .'3x/3)  R30'  sunerst  rueture,  and  oxidation  at  room 
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20.  Abstract  (continued) 


•  temperature  saturates  at  a  stage  which  incorporates  V)  monolayers  of  nickel  in 
a  stoichiometric  amorphous  film  of  NiO  whereas  at  temperatures  higher  than 
^200JC  thicker  oxide  films  are  produced.  Our  study  indicates  that  each  oxygen 
atom  adsorbed  on  the  N'i(ill)  surface  interacts  with  and  relaxes  three  nearest 
neighbor  nickel  atoms,  and  after  saturation  of  the  relaxation,  oxidation  of 
three  monolayers  takes  place  abruptly  after  which  the  oxide  layer  on  the  sur¬ 
face  apparently  blocks  further  reaction. 
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Abstract.  Surface  atomic  structures  of  clean,  oxygen-adsorbed,  and  oxidized  (111) 
nickel  have  been  studied  quantitatively  by  using  MeV  ion  scattering  in  combination 
with  AES  and  MEED.  We  show  that;  the  clean  (111)  nickel  surface  has  the  bulk-like 
structure  with  reconstruction  or  relaxation  less  than  0.02  X,  the  surface  thermal 

vibration  amplitude  is  enhanced  by  ^20  %  compared  to  the  bulk  value,  adsorbed 

o 

oxygen  results  in  surface  lattice  expansion  by  “vO.15  A  which  is  closely  correlated 

o 

to  the  p(2x2)  or  (/3x/3)R30  superstructure,  and  oxidation  at  room  temperature 
saturates  at  a  stage  which  incorporates  ^3  monolayers  of  nickel  in  a  stoichiometric 
amorphous  film  of  NiO  whereas  at  temperatures  higher  than  *\/2OO0C  thicker  oxide  films 
are  produced.  Our  study  indicates  that  each  oxygen  atom  adsorbed  on  the  Ni(lll) 
surface  interacts  with  and  relaxes  three  nearest  neighbor  nickel  atoms,  and  after 
saturation  of  the  relaxation,  oxidation  of  three  monolayers  takes  place  abruptly 
after  which  the  oxide  layer  on  the  surface  apparently  blocks  further  reaction. 
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1.  Introduction. 

The  field  of  surface  crystallography  is  still  in  its  infancy.  Accurate  atomic 
positions  are  often  not  known  even  on  rather  simple-structured  surfaces  which  show 
a  bulk-like  (lxl)  diffraction  pattern.  The  agreement  between  measurements  and 
calculations  of  LEED  I-V  curves,  which  has  been  a  major  tool  for  surface  structure 
investigation,  is  often  unsatisfactory*  in  uniquely  defining  the  surface  structure. 
The  situation  is  even  worse  when  the  symmetry  of  the  structure  is  lowered  by 
reconstruction  or  adsorption- induced  superstructures.  Typically,  for  ordered  species 
the  adsorbate  binding  site  and  its  height  over  the  substrate  are  used  as  model 
parameters  for  analysis  of  model  structures  using  LEED  data  usually  under  an  assump¬ 
tion  of  unstrained  substrate  surface  structure.  The  success  of  this  approach  is 
often  limited  by  the  number  of  model  structures  that  can  be  analyzed.  On  the  other 
hand,  there  is  growing  evidence  that  the  geometry  of  the  substrate  may  be  changed 
appreciably  due  to  gas  adsorption.  In  the  particular  case  of  Ni(lll)  surface, 
recent  LEED  measurements  of  temperature  vs.  coverage  phase  diagram  for  the  0/Ni(lll) 
system^  and  high  resolution  electron  energy  loss  measurements  for  the  H/Ni(lll) 
system  have  both  suggested  that  the  Ni(lll)  surface  lattice  might  be  reconstructed 

due  to  adsorption.  We  have  presented  definitive  evidence  of  surface  relaxation  in 

3) 

the  0/Ni(lil)  system  by  high  energy  ion  scattering  measurements  . 

In  the  present  paper,  we  describe  detailed  experimental  results  of  KeV  He  ion 
scattering  for  the  0/Ni(lll)  system,  compare  them  to  the  results  of  computer  simula¬ 
tion,  and  discuss  possible  structure  models  for  the  oxygen-induced  p(2x2)  and 
(/5x/3)R30#  superstructures  and  also  the  interface  structure  between  the  crystalline 
nickel  substrate  and  the  oxide  overlayer.  The  usefulness  of  the  so-called  surface 
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peak  intensity,  vhich  is  measured  in  an  aligned  backscattering  arrangement,  has 

k) 

been  well  demonstrated  elsewhere  as  a  structure-sensitive  probe  for  surface 
studies.  In  addition,  the  ion  scattering  measurements  combined  with  in-situ  AES 
allows  absolute  oxygen  coverage  determinations  over  a  wide  range  from  submonolayer 
to  over  10  monolayers  coverage. 

2.  Experimental. 

Clean  Ni(lll)  surfaces  were  prepared  in  UKV  by  oxidation  and  reduction  cycles 
and  argon  ion  bombardment  followed  by  annealing.  Prolonged  slow  argon  ion  bombard¬ 
ment  with  the  target  at  ''■'600° C  was  especially  effective  to  remove  the  impurities 
carbon  and  sulfur  vhich  are  common  segregations  on  metal  surfaces. 
ar.w ...Tgc*ii  aan-aauicT arjunr  waKt *****  mmmmmmmm*.  The  cleanliness  and  symmetry  of 

•the  target  surface  were  examined  by  in-situ  AES  and  5  keV  electron  diffraction 
(MEED;.  After  the  cleaning  proctdure,  we  measured  a  typical  clean  nickel  Auger 
spectrum  with  no  impurities  exceeding  0.1  %  in  the  surface  atomic  composition.  The 
oxygen  coverage  in  the  initial  stage  of  adsorption  (  ^1  monolayer  =  1.86x10 
atoms/cm  )  was  monitored  by  AES,  whose  absolute  sensitivity  was  calibrated  by  ion 
backscattering  measurements,  and  greater  oxygen  coverages  (  ^1  monolayer)  in  the 
later  adsorption  or  oxidation  stage  were  measured  directly  by  backscattering.  In 
accordance  with  previous  studies^’^,  ue  found  two  well  developed  superstructure 
diffraction  patterns,  p(2x2)  and  (/3x*/5)P3C°  ,  in  the  course  of  oxygen  adsorption 
at  room  temperature.  The  diffraction  pattern  gradually  disappeared  after  the  oxygen 
coverage  of  ^1/3  monolayers  indicating  the  onset  of  amorphous  oxide  film  formation. 

The  He+  ion  beam  of  0.?  to  2  MeV  produced  by  SUKY-Dynanitron  was  introduced 
to  the  UKV  target  chamber  through  two  1  nm  orifices.  Scattered  particles  were 
analyzed  by  a  Si  detector  of  surface-barrier  type,  with  sn  energy  resolution  of 
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^20  keV  located  in  a  grazing  exit-angle  geometry  (^5°  to  the  surface)  to 
increase  the  surface  sensitivity.  The  target  orientation  with  respect  to  the 
incident  beam  vas  aligned  in  either  the  normal  or  non-nornal  axial  channeling 
direction  (see  a  marble  model  in  Fig.l).  Due  to  shadowing  of  the  underlying 
atoms  in  a  string  by  the  surface  atom,  we  can  clearly  isolate  and  evaluate  the 


surface  peak  (SP)  in  the  backscattering  spectrum  as  shown  in  Fig. 2.  The  oxygen 
peak  intensity  vas  evaluated  after  subtracting  the  Hi  background,  which  was  estima  ed 
by  taking  the  average  of  the  counts  in  each  10  channels  on  both  sides  of  the  oxygen 


Thus,  we  could  determine  the  absolute  oxygen  coverage  down  to  tl  nonolayer  with  an 


accuracy  of  +15  %.  The  areal  intensity  of  the  S?  was  evaluated  by  using  a  standard 

background  subtraction  method^  and  compared  to  the  random  scattering  yield  to 
determine  the  number  of  backscattering  Ki  atoms  (atoms/row). 

In  measurements  of  the  angular  distribution  of  the  SP  intensity,  we  obtained 

a  series  of  spectra  as  shown  representatively  in  Fig. 3(a).  We  define  the  bulk 

yield  and  the  surface  yield  as  the  integrated  counts  in  the  segments  labelled 

"B"  and  "S",  respectively,  in  the  spectrum.  If  we  plot  these  yields  as  a  function 

of  the  incident  angle  we  observe  the  familiar  channeling  dips  as  shown  in  Fig. 3(b). 

The  bulk  yield  shows  a  symmetric  shape  and  its  centroid  can  be  used  to  locate  the 

exact  channeling  direction  (a  <110>  axis  in  the  case  of  Fig. 3).  In  contrast, 

the  surface  yield  shows  some  asymmetric  distribution  which  will  be  discussed  in 

detail  and  related  to  the  surface  atomic  structure.  As  shown  in  Fig. 3(a),  it 

gradually  becomes  harder  to  isolate  the  SP  in  the  spectrum  with  increasing 

misalignment.  Therefore,  detailed  spectral  intensities  were  compared  with 

calculation  only  within  an  angular  range  of  'til  deg.  from  the  axial  channeling 

direction.  We  can  still  isolate  the  SP  in  this  range,  and  we  applied  a  similar 
background  subtraction  method  as  that  described  in  ref. 6  though  this  method  becomes 
increasingly  inaccurate  with  increasing  minimum  yield  behind  the  SP  as  shown  in 
spectra  Q  and  Q  in  Fig. 3(a).  The  fluence  of  He+  ion  beam  vas  M.  yC/mn^  (n.6xl0 
ions/cm2)  per  each  spectral  measurement,  which  we  determined  to  be  far  (^2  orders) 


I 


W— far  fceiow  the  critical  value  at  vhich  the  effect  of  radiation  damage  or 

beam  induced  desorption  becomes  significant.  The  background  pressure  in  the  target 

chamber  did  not  exceed  5x10”"^  Torr  during  measurements. 


3.  Computer  simulation. 

The  SP  intensity  of  nickel  was  calculated  by  the  nuclear  encounter 

7) 

probability  approach  assuming  a  number  of  surface  atomic  structure  and  dynamics 

5  7-9' 

models.  Details  of  the  calculation  method  have  been  discussed  in  the  literature  ’ 

In  brief,  the  paths  of  charged  particles  vith  given  energy  and  direction,  vhich  impinge 

uniformly  relative  to  a  single  string  of  target  atoms,  vere  traced  in  a  stepvise 

manner  by  using  a  binary  collision  approximation.  The  consequent  nuclear  encounter 

probability  for  each  atom  in  the  string  was  summed  to  obtain  the  SP  intensity. 

The  one-dimensional  rms  amplitude  for  bulk  nickel  atoms  vas  calculated  to  be 

0.068  i  assuming  the  target  temperature  of  300  ^310  K  during  measurements  and  a 

bulk  Debye  temperature  of  U20  K.  This  value  was  used  throughout  the  present  study 

for  the  bulk  amplitude,  and  the  surface  enhancement  of  thermal  vibration  vas 

introduced  through  an  exponential  decay  model  similar  to  the  calculation  by 

Stensgaard  et  al^,  i.e.  the  enhancement  vas  assumed  to  decay  exponentially  into 

o 

the  bulk  vith  a  decay  constant  of  3.52  A  (lattice  constant  of  nickel).  Possible 
anisotropy  in  the  surface  thermal  vibration  vas  not  taken  into  account.  We  used 
a  coefficient  of  0.2  for  the  correlation  of  thermal  vibration  though  an  exact 
value  for  nickel  is  not  knovn.  Based  on  experimental  phonon  frequency  distributions, 
Jackson  et  al10’11^  have  calculated  the  correlation  coefficients  for  Mo  and  Ifb  to 


l 


a 


be  between  ''-0.2  and  ^.3  at  ^300  K  depending  on  the  interatomic  distance.  These 
values  are  slightly  lover  than  the  Debye  model  prediction.  We  calculated  the 
dependence  of  the  SP  intensity  on  the, correlation  coefficient  and  confirmed  that 
the  calculated  result  does  not  differ  more  than  ±5  %  if  we  choose  a  value  between 
0.1  and  0.3. 


Shadowing  of  underlying  nickel  atoms  by  adsorbed  oxygen  was  estimated  using  a 
similar  method  as  that  described  in  ref. 9-  The  effect  is  maximized  when  the  adsorbed 
oxygen  sits  exactly  on  top  of  the  Si  string  and  decreases  steeply  with  the  misfit 
distance  between  the  oxygen  site  and  the  string  axis.  In  order  to  estimate  the 
“maximum  effect,  we  assumed  the  3-fold  oxygen  adsorption  sites  (A  and  B  in  Fig.l). 


t  ,  .ft 

The  0-Ni  distance  was  assumed  to  be  (2.03+l*20)  A  for  site-A  and  (2x2.03+1.20)  A 

for  site-B,  where  the  distance  of  oxygen  atom  from  the  surface  (1.20  £.)  was  that 

12 ) 

obtained  from  LEED  analysis  .  The  results  of  simulation,  however,  is  not  very 

q  \  +  . 

sensitive  to  this  0-Ni  distance^  .  The  single  string  calculation  for  1  MeV  Re  ions 

parallel  to  the  <111>  direction  resulted  in  the  shadowing  effects  of  ^0.17  atoms/row 

and  M).20  atoms/row  for  one  monolayer  of  site-A  and  site-B  oxygen,  respectively.  Since 

osar  the  oxygen  coverage  for  the  p(2x2)  or  (/3x/3)R30  superstructure  does  not 

exceed  1/3  monolayers  as  will  be  discussed  later,  the  net  decrease  of  the  SP 

intensity  by  the  shadowing  effect  must  be  averaged  and  is  thereforeM.^  or  less  of  the 

single  string  calculation  result  ,  i.e.  the  shadowing  effect  by  oxygen  in  the 

superstructures  turns  out  to  be  '''0.07  atoms /row  at  most.  This  is  small  enough 

to  be  neglected  compared  to  the  systematic  error  of  measurements  which  we  estimated 

to  be  +5  %.  If  ve  assume  the  position  of  oxygen  atom  as  described  above,  no 

shadowing  effect  is  expected  in  the  <110>  incidence  case  because  the  oxygen  atom 

is  away  from  the  <110>  axis  by  'u0.U8  %.  We  have  earlier  demonstarted  that  the 

9) 

shadowing  effect  from  an  amorphous  overlayer  is  negligible  .  The  conclusion  is 
that  ve  can  safely  neglect  the  shadowing  effect  by  adsorbed  oxygen  or  amorphous 
oxide  film  in  the  following  discussion. 
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In  simulations  of  the  angular  scan  of  the  SP  intensity  for  oxygen  adsorbed 

surfaces  we  have  also  used  a  single  string  approximation.  This  becomes  inaccurate 

as  the  beam  direction  deviates  from  the  channeling  directions.  We  have  previously 

9) 

found  by  using  a  multi-string  simulation  method  that  small  contributions  from 
deeper  atoms  begin  to  be  significant  when  the  incident  angle  deviates  more  than 
about  1  deg.  from  the  channeling  direction.  However*  for  deviation  less  than 
±1  deg.,  the  contribution  from  all  atoms  deeper  than  the  12th  atom  in  a  string 
sums  up  to  less  than  2  %  of  the  total  yield.  Thus  we  can  neglect  the  error  due  to  the 
single  string  model  within  the  angular  range  of  ±1  deg.  from  the  channeling 
direction. 


^4.  Results  and  Discussion. 

U.l  Clean  surface. 

The  clean  (ixl)  surface  shows  a  dependence  of  the  SP  intensity  on  incident 

He+  ion  energy  as  shown  in  Fig.U.  Since  our  technique  is  sensitive  to  atomic 

displacements  perpendicular  to  the  incident  direction,  normal  incidence  measurements 

(a)  give  information  about  lateral,  or  in-plane,  displacements  of  surface  atoms, 

and  non-normal  incidence  measurements  (b)  are  sensitive  to  both  lateral  and 

longitudinal,  or  out  of  plane,  displacements.  The  monotonic  increase  of  the  SP 

intensity  in  both  figures  corresponds  to  the  decrease  of  shadowing  effect  with 

-1/2 

increasing  incident  energy  and  follows  nicely  the  expected  E  dependence  of 
the  so-called  shadow  cone  radius.  Solid  curves  (a)  and  (b)  in  both  figures  are 
results  of  computer  simulation.  Both  curves  assume  a  bulk-like  atomic  structure. 
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and  curve  (a)  assumes  a  bulk-like  thermal  vibration  while  curve  (b)  assures  a 
small  enhancement  (  n. 20  %)  in  the  surface  thermal  vibration  amplitude.  The  agree¬ 
ment  of  curve  (b)  with  the  experimental  plots  in  both  figures,  together  with  the 
(lxl )  diffraction  pattern,  strongly  suggests  that  the  clean  (ill)  nickel  surface 
has  the  bulk  termination  structure.  This  agreement  also  suggests  that  anisotropy 
in  the  thermal  vibration  is  small.  In  fact  simulations  indicate  that  anisotropy 

in  the  ms  amplitude  by  as  small  as  0.00U  A  should  appear  as  a  detectable 
variation  in  the  S?  intensity  in  the  two  different  incident  directions.  This  result 
is  in  reasonable  agreement  with  Allen  and  de  Vette!-2a',  who  calculated  the  transverse 
and  longitudinal  enhancements  in  the  mean  square  amplitude  to  be  a  factor  of  1.30 
and  1.86,  respectively,  in  the  high  temperature  limit  approximation  assuming  force 
constants  at  the  surface  equal  to  those  in  the  bulk.  Enhancement  of  the  surface 
thermal  vibration  assumed  in  curve  (b)  corresponds  to  a  surface  Debye  temperature,  6jj, 
of  about  330  K,  which  is  in  a  good  agreement  with  the  DEED  study!-2)  though  another 
T/FFD  study  by  Unertl  cited  in  ref. 9a  suggests  a  lower  value  of  ^150  K. 

The  symmetric  angular  distribution  of  the  <110>  SP  intensity  shown  in  Fig. 5 

confirms  the  above  conclusion.  Solid  and  dashed  curves  in  this  figure  shows  the 

results  of  computer  simulation  in  which  the  bulk-like  structure  and  an  expansion 

of  the  first  monolayer  by  0.05  X  are  assumed,  respectively.  The  same  enhanced 

surface  thermal  vibration  as  in  Fig.!*  is  also  assumed.  The  simulation  result  shows 

a  significant  asymmetry  in  the  angular  distribution  even  for  relaxation  of  the 

first  monolayer  of  only  0.02  A.  Consequently  the  observed  symmetric  distribution 

confirms  that  the  surface  structure  is  bulk-like.  A  slightly  narrower  distribution 
of  the  experimental  plots  compared  to  the  solid  curve  suggests  a  possible  higher 
vibrational  amplitude  hence  lower  value  of  than  that  used.  Adoption  of  a  value 
that  fits  the  apparent  width  however  leads  to  a  too  high  value  of  the  S?  intensity 
in  the  aligned  direction.  We  feel  that  the  apparently  narrower  experimental  distri¬ 
bution  is  a  result  of  the  larger  uncertainty  and  probable  insufficient  magnitude  of 
the  background  subtraction  at  the  larger  incidence  angles. 

However,  there  is  no  doubt  about  the  symmetry  of  the  angular  distribution.  We 
can  conclude  from  this,  from  the  agreement  between  the  calculated  and  measured 
SP  intensity  in  the  vicinity  of  the  <110>  axis  and  from  the  energy  dependence 
shown  in  Fig.U  that  the  clean  (ill)  nickel  surface  has  a  bulk-like  atomic  structure 
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vith  relaxation  or  reconstruction  less  than  0.02  end  vith  ^20  /»  nearly  isotropic 

enhancement  of  surface  thermal  vibration  amplitude.  This  results  confirms  the 

j.2  13 ) 

conclusion  by  LEZD  studies^  ’  in  a  straightforward  end  somewhat  more 
quantitative  manner. 


k.2  Oxygen  adsorption. 

The  S?  intensities  for  <111>  and  <11Q>  beam  directions  as  a  function  of  exposure 
of  the  crystal  surface  to  0 g  gas  at  room  temperature  (30015  K)  are  shown  in  rig. 6  . 

A  striking  difference  between  the  two  incident  directions  at  low  oxygen  exposure 
indicates  relaxation  of  surface  atoms  in  a  direction  perpendicular  to  the  surface. 
Also  clear  steps  corresponding  to  superstructure  MEED  patterns  as  shown  in  the 

figure  suggest  that  the  relaxation  is  closely  related  to  the  superstructure.  The 

absolute  oxygen  coverage  measured  at  the  exposures  indicated  by  arrows  in  Fig. 6  were 
0.23±.0v  monolayers  for  the 

ip  -i  p(2x2)  structure  and  0.31±.05  monolayers  for  the  (/3x/3)R30  structure. 

The  oxygen  coverage  increases  to  "MD. 36  monolayers  atVLO  Langmuir  exposure  where  an 
abrupt  increase  in  the  S?  intensity  occurs  in  both  incident  directions.  This 
increase  is  accompanied  vith  rapid  disappearance  of  the  MEED  pattern,  indicating 
the  onset  of  amorphous  oxide  film  formation.  These  results  are  consistent  vith 
LEED-AES  studies  5  .  The  results  shown  in  Fig. 6  constitute  a  direct  observation 

of  surface  relaxation  of  nickel  substrate  atoms  connected  with  formation  of  the 
superstructures.  However,  since  Fig. 6  alone  cannot  lead  us  to  a  unique  structure 
model,  we  analysed  the  energy  and  angular  dependences  of  the  SP  intensity  in 


detail. 


Figure  7  shows  the  energy  dependence  of  the  <110>  S?  intensity  at  the  stage 
of  the  (/3x/3)R30C  structure.  Solid  curves  were  calculated  assuming  uniform 
relaxation  cf  the  first  monolayer  in  the  perpendicular  direction  to  the  surface 
by  different  amounts  varying  from  0.0  &  for  the  bottom  curve  to  0.52  %  for  the  top 
curve.  Since  the  <110>  direction  is  35.26  deg.  from  the  normal  direction,  these 
relaxations  can  be  seen  by  the  incident  beam  as  displacements  of  the  first  atomic 
layer  from  the  bulk  lattice  position  from  0.0  X  to  0.30  A  in  steps  of  0.05  The 

observed  SP  intensity  lies  close  to  or  slightly  below  the  third  curve  from  the 

o 

bottom.  This  suggests  that  the  relaxation  of  the  first  monolayer  is  0.151.02  A. 

In  this  measurement  outward  and  inward  relaxation  play  exactly  the  same  role,  so 
ve  cannot  distinguish  between  expansion  and  contraction  of  the  surface.  Also  the 
relaxation  may  not  be  restricted  to  the  first  monolayer  since  some  combination 
of  surface  and/or  subsurface  strains  would  yield  almost  the  sane  simulation 
result  as  shown  in  Fig. 7- 

Figure  8(a)  and  (b)  compare  the  results  of  computer  simulations  with  experi¬ 
mental  angular  scans  of  the  <110>  SP  intensity  for  the  p(2x2)  and  (/3x»^3)n3O0 
structures,  respectively.  In  the  calculations  we  assumed  uniform  expansion  of  the 
first  monolayer  by  various  amounts  as  indicated  for  each  curve  in  the  figure.  If 
contraction  instead  cf  expansion  is  assumed,  the  asymmetry  and  shift  of  the  minimum 
point  are  reversed.  This  is  not  consistent  with  observation,  so  we  can  rule  out  surf 
contraction  for  both  superstructures.  In  addition,  since  assumptions  of  relaxation 
in  the  second  or  subsurface  layers  do  not  yield  proper  angular  distributions ,  we  can 
essentially  limit  the  relaxed  atoms  to  the  first  layer. 

The  experimental  plots  in  Fig. 8(a)  apparently  agree  with  an  expansion  of  M).12 
On  the  other  hand,  the  plots  in  Fig. 8(b)  are  much  closer  to  the  simulation  result 
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vhich  assumes  a  uniform  expansion  of  the  first  layer  by  0.15  A  in  a  good  agreement 

with  the  energy  dependence  shown  in  Fig. 7.  In  order  to  further  investigate  the 

situation  for  the  p(2x2)  structure,  we  calculated  some  non-uniform  expansion  models 

as  shown  in  Fig. 9.  The  best  fit  is  with  curve  (b)  which  assumes  that  3/1*  of  the 

nickel  atoms  in  the  first  monolayer  are  expanded  by  0.15  &  ana  1/U  are  not  expanded 

3) 

at  all.  Similar  non-uniform  expansion  models  were  also  compared  to  the  experimental 
plots  for  the  (/3x/3)R30°  structure,  however,  in  no  case  did  we  obtain  better  agree¬ 
ment  than  shown  for  the  solid  curve  in  Fig. 8(b). 

In  the  angular  scan  computer  simulations  we  used  the  same  enhanced  thermal 
vibration  model  as  determined  in  li.l  for  the  clean  surface.  This  approach  is 
indicated  by  the  constancy  of  the  observed  S?  intensity  with  oxygen  coverage 

up  to  '■'10  Langmuir  for  the  <111>  beam  direeiion(Fig.6 ) .  Significant  decrease  in  the 
thermal  vibrational  amplitude  would  result  in  decrease  (or  in  the  unlikely  event  of 
a  vibration  increase,  in  increase)  in  this  direction,  independent  of  relaxation  of 
the  surface  layer. 

All  of  our  results  taken  together;  that  the  oxygen  coverage  is  Vl/k  monolayers 
for  the  p(2x2)  and  'tl/3  monolayers  for  the  (/3x/3;E30°  structures,  that  3A  cf  the 
first  monolayer  is  relaxed  for  the  p(2x2)  ar.d  all  cf  the  first  monolayer  is  relaxed 
for  the  (/3x/3)R30°  structure  and  the  LEED  evidence  that  the  3-fold  symmetry  site 
is  most  likely  for  oxygen  adsorption  on  (ill)  nickel  surface  leads  to  the  atomic 
structure  models  shown  in  Fig. 10  for  the  (lxl),  p(2x2)  and  (/Jx*/3)R30°  structures, 
Oxygen  atoms  adsorbed  on  3-fold  sites  (site-B  in  Fig.l)  are  shown  by  large  shadowed 
circles,  and  small  closed  and  open  circles  represent  the  expanded  and  non-expanded 
(bulk-like)  nickel  atoms  in  the  first  monolayer,  respectively.  Neither  the  LEED  nor 
ion  scattering  results  differentiate  between  the  A  or  B  adsorption  sites.  For  both 
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structures  only  one  type  of  site  is  occupied-perhaps  electronic  structure 
calculations  can  be  used  to  identify  the  lowest  energy  of  the  tvo  possible 
configurations.  The  transition  from  the  p(2x2)  to  the  (/3x/3)R30°  structure  in 
this  model  involves  motion  of  1/2  of  the  oxygen  atoms  to  other  bonding  sites. 
Formation  of  the  p(2x2)  structure  over  the  entire  surface  before  the  phase  change 
to  the  (/5x/3)R23°  structure  implies  that  it  is  a  lower  energy  configuration  than 
islands  of  the  (/3x/3)R30°  structure.  In  any  case  our  atomic  model  suggests  the 
Ni^O  cluster  as  an  electronically  stable  composition  on  the  oxygen  adsorbed  (ill) 
nickel  surface. 

1».3  Oxidation. 

When  the  crystal  surface  is  exposed  to  0^  gas  at  room  temperature  both  of  the 

<111>  and  <110>  S?  intensities  increase  abruptly  at  an  exposure  of  'VIO  Langmuir, 

where  the  oxygen  coverage  was  measured  to  be  ^.36  monolayers.  The  increase  in  the 

SP  intensities  apparently  saturates  after  an  exposure  of  about  100  Langmuir  as  shown 

in  Fig. 6.  The  saturation  increase  is  M.  atoms/row  for  the  <111>  S?  intensity  and 

n»3  atoms /row  for  the  <110>  S?  intensity.  Since  the  incident  beam  sees  directly  top 

three  monolayers  in  a  <U1>  direction  and  only  the  first  monolayer  in  the  <110> 

incident  direction,  both  of  these  saturation  increases  correspond  to  release  cf  the 

first  three  monolayers  of  nickel  atoms  from  registration  with  the  lattice.  In  detail 

this  result  shows  that  the  first  three  monolayers  are  displaced  an  amount  large 

compared  to  the  shadow  cone  radius  at  2  MeV  (  'vo.Ofi  ^  in  the  <110>  incidence  and 
M).l  A  in  the  <111>  incidence)  and  are  not  correlated  vith  each  other.  This,  together 

with  the  disappearance  of ' the  MEED  pattern  at  this  stage,  suggests  an  amorphous  oxide 

layer  formation.  Measurements  of  the  angular  distribution  of  the  nickel  and  oxygen 

scattering  intensity  support  'this  conclusion.  Figure  11  shows  the  energy  dependenc 

of  the  SP  intensity  for  the  oxygen  saturated  surface  at 


room  temperature.  We  can  derive  some  information  about  the  interface  atomic 

structure  form  these  measurements  by  using  the  concept  of  the  shadow  cone  radius. 

Figure  11(b)  shows  that  the  difference  between  the  "clean"  and  "oxygen  saturated" 

SP  intensities  is  ^3  atons/row  for  <110>  incidence  over  the  whole  energy  range 

studied,  which  confirms  large  displacements  of  the  first  three  monolayers.  On  the 

other  hand,  Fig. 11(a)  shows  that  the  similar  difference  in  the  <111>  incidence 

case  is  atom/row  in  the  high  energy  range  but  decreases  to  t«0.7  atoms/row  in 

the  low  energy  range.  The  shadow  cone  radii  at  0.6  MeV  at  the  position  of  the 

second  atom  in  each  row  are  ^0.116  A  for  <11Q>  incidence  and  ru0.l8l  %  for  <111> 

incidence.  Therefore,  Fig. 11(a)  and  (b)  indicate  that  displacements  of  the  nickel 

o 

in  the  top  three  layers  are  veil  above  'tO.12  A  in  a  direction  perpendicular  to 

the  <110>  direction.  However,  the  displacements  of  some  of  the  atoms  in  the  three 

layers,  probably  the  layer  nearest  the  interface,  are  %0.l8  A  vhen  viewed  in  the  direc 
normal 

-  to  the  surface.  This  suggests  that  the  oxide  overlayer  is  not  really  amorphous 
throughout  but  some  registration  still  exists  to  the  underlying  crystal,  whose 
structure  is  not  distorted  within  the  accuracy  of  the  measurement  ('VO.O?  £) . 


The  absolute  oxygen  coverage  measurement  at  the  saturation  stage  resulted 


in  2.8±.t*  monolayers  confirming  formation  of  n-3  monolayers  of  oxide  film  giving 

a  stoichiometry  of  "iO.  This  agrees  fairly  veil  with  a  nuclear  microanalysis 
measurements  by  Norton  et  ai ^ ' ,  however,  it  disagrees  with  LEED-AZS  measurement 

which  concluded  the  saturation  coverage  to  be  1.6  layers  of  KiO  at  room 


temperature.  We  believe  this  disagreement  is  due  to  the  difficulty  of  AES  in 


performing  absolute  coverage  measurements,  especially  in  a  range  thicker  than 
1  monolayer.  Previous  studies  cf  the  oxidation  of  the  nickel  surface  using  LEEIA- 
AES^^  and  ion  scattering^ ^  have  concluded  that  the  saturation  thickness  is  \2 


monolayers  for  both  the  (100)  and  (110)  surfaces,  respectively.  If  this  is  true. 
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it  is  interesting  to  note  that  both  of  these  surfaces  expose  the  top  two  layers 
directly  to  the  vacuum. 


l*.l*  High  temperature  measurements . 

When  the  crystal  was  kept  at  temperatures  higher  than  ^350  K  and  less  than 
*'J*7o  K  during  the  oxygen  exposure  and  measurements,  we  could  discern  a  weak 
diffraction  pattern  corresponding  to  the  p(2x2)  structure  in  the  initial  stage 
of  oxygen  adsorption,  but  the  (/3x/3)R30°  structure  was  not  detected.  Neither  of 
these  structures  were  found  if  the  crystal  temperature  exceeded  n«U70  K.  Instead, 
we  found  a  diffuse  and  streaky  diffraction  feature  after  an  exposure  of  MO  L, 
which  seemed  to  become  sharper  and  stronger  with  increased  oxygen  exposure  and 

increased  crystal  temperature  which  indicates  crystali2ation  of  the  KiO  film. 

1  1& ) 

These  KEED  results  are  in  a  good  agreement  with  previous  LEED  studies  * 

Variation  of  the  <110>  SP  intensity  is  shown  in  Fig. 12  for  1*20  K  and  520  E 
measurements.  The  SP  intensity  for  2  MeV  ions  for  the  clean  surface  at  these 
temperatures  were  U.32  and  U.i*8  atoms/rov,  respectively.  The  difference  between 
these  values  and  that  at  room  temperature  shown  in  Fig.U(b)  is  due  to  the  increase 
in  the  thermal  vibration  and  agrees  well  with  computer  simulations.  In  case  of 
1*20  K  measurements,  the  initial  increase  in  the  SP  intensity  in  the  exposure 
region  lower  than  MO  L  is  about  the  same  as  that  corresponding  to  the  p(2x2) 
structure  in  Fig. 6.  Although  detailed  angular  or  energy  distribution  measurements 
have  not  been  carried  out,  the  surface  strucutre  is  probably  similar  to  that 
shown  in  Fig. 10.  As  observed  for  lover  temperature  measurements  the  SP  intensity 
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increases  at  aroud  10  L.  However,  the  oxygen  coverage  at  10  L  was  measured  at 
^0.31  monolayers  (vs.  'VD.36  monolayers  at  room  temperature)  and  the  rate  of  increase 
in  the  SP  intensity  is  lover  than  at  room  temperature.  We  believe  these  differences 
arise  from  a  lover  sticking  probability  of  oxygen  on  the  (ill)  surface  at  the 
high  temperature.  The  apparent  saturation  level  at  1+20  K  is  M3  atoas/rov  which 
agrees  with  room  temperature  measurements.  At  520  K,  the  initial  increase  in  the 
SP  intensity  is  slightly  higher  than  the  room  temperature  or  1+20  K  measurements. 

This,  together  with  no  diffraction  feature  observed  at  this  temperature,  is 
consistent  with  the  suggestion  from  LEED  and  ion  scattering  studies"'"** that 
adsorbed  oxygen  reacts  with  nickel  and  makes  small  clusters  of  NiO.  The  possible 
intermediate  decrease  in  the  S?  intensity  at  ^8  L  may  be  explained  by  agglomeration 
of  small  clusters  to  larger  islands.  Increase  of  the  SP  intensity  after  MOL 
exposure  is  also  slow  in  this  case  but  continues  without  saturation.  We  have 
measured  up  to  an  exposure  of  900  L  where  the  increase  of  SP  intensity  was  Ml6 
atoms/row,  which  indicates  oxide  film  formation  involving  more  than  M6  monolayers 
of  nickel.  The  rate  of  increase  shows  a  tendency  to  decrease  with  increasing  thickness 
suggesting  that  the  oxidation  may  be  limited  by  diffusion.  Although  comparison 
of  aligned  and  random  spectra  in  this  stage  showed  a  significant  difference  in 
the  nickel  and  oxygen  peak  intensity  indicating  crystalline  structure  of  the 
oxide  film,  the  ratio  of  the  aligned  peak  intensity  to  the  random  (  %)  is 

evidence  of  considerable  disorder  in  the  film.  The  temperature  dependence  of  the 
oxidation  presented  here  is  in  a  good  agreement  with  previous  studies. 


5.  Summary. 


We  have  applied  high  energy  ion  scattering  to  a  study  of  clean,  oxygen- 
adsorbed,  and  oxidized  (ill)  surface  of  nickel  and  shown  that; 

1)  The  clean  surface  has  the  bulk-like  structure  with  relaxation  or  reconstruction 
less  than  ^0.02  jl. 

2)  The  surface  thermal  vibration  of  the  clean  surface  is  enhanced  by  %20  % 
compared  to  the  bulk, 

3)  The  p(2x2)  structure  induced  by  adsorption  of  M.A  monolayers  of  oxygen  is 

accompanied  by  outward  relaxation  (lattice  expansion-)  of  *u3A  of  the  nickel 

o 

atoms  in  the  first  monclayerby  ''-O.IJ  A  perpendicular  to  the  surface. 

U)  The  (/3x/3)R30°  structure  induced  by  adsorption  of  n, 1/3  monolayers  of  oxygen 

is  accompanied  by  outward  relaxation  of  all  of  the  atoms  in  the  first  monolayer 

0 

of  the  nickel  crystal  by  'tO.15  A. 

5)  Oxidation  of  the  (ill)  nickel  surface  takes  place  at  '^0.36  monolayers  coverage 
at  room  temperature  and  at  slightly  lover  coverage  at  higher  temperatures. 

6)  Saturation  of  the  oxide  film  thickness  involves  ^3  monolayers  of  nickel 
at  temperatures  lower  than  'A70  K.  At  higher  temperatures  there  is  r.o 
saturation  in  the  exposure  range  studied. 

7)  The  oxide  composition  is  KiO.  It  is  amorphous  with  some  indication  of 
registration  at  the  interface  with  the  substrate  crystal  at  room  temperature , 
and  crystalline  at  higher  temperatures. 

8)  The  interface  structure  is  very  sharp  with  distortion  of  the  substrate  lattice 
less  than  'V0.05  X  at  room  temperature. 
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These  results  are  consistent  with  previous  studies  ggssssSasaa 5 — -•*  1 — ■ 

n  ■■  2^3  provide  complementary 

new  information  on  the  surface  structure.  A  detailed  interface  study  of  the 
epitaxial  system  such  as  JJiO/Ni ( 111 )  presented  here  briefly  is  an  interesting 
topic  for  high  energy  ion  scattering,  however,  it  is  beyond  the  scope  of  this 
paper. 

Acknowledgement  is  due  to  L.C. Feldman,  Eell  Telephone  Laboratories, 
and  T.Gustaffscn  and  E.W. Plummer,  University  of  Pennsylvania,  for  invaluable 
suggestions  and  stimulating  discussions.  This  work  was  supported  by  the  Office 
of  Kaval  Research  under  Contract  No .IJ0001^-78-C-06l6.  One  of  us  (ET)  acknowledges 
the  support  from  the  Office  of  Naval  Research  under  Contract  No .N0001U-79-C-0991. 


-18- 


References. 

1)  L.D.Roeiofs,  A.R.Kortan,  T.L. Einstein  and  R.L.Park; 

in  J. Vac. Sci .Technol.  18  (1981 )  **92. 

2)  W.Ho,  K. J.DiNardo  and  E.W-.Plumer ;  J.Vac. Sci. Technol.  17.  (l93o)  130 

3)  T.llarusawa,  W.M.  Gibson  and  E.Tornqvist;  Phys.  Rev.  Letters  1*7.  (1981)  1*17  • 
1*)  E. Bogh  in  Channeling:  Theory,  Observation  and  Applications ,  ed.  by 

D.V.Korgan,  John  Wiley  &  Sons,  Hew  York  (1973). 

L.C.Feldran  in  Surface  Science:  Recent  Progress  and  Perspectives,  ed.  by 
R.Vanselov,  CRC  Press,  Ohio  (1980). 

5)  H.Tbach  and  D.Bruchrann;  Phys.  Rev.  Letters  U*^  (1980)  36. 

6)  I.Stensgaard,  L.C.Feldran  and  P.J.Silveman;  Surface  Sci.  77  (1978)  513- 

7)  J.H. Barrett ;  Phys. Rev.  B3  (1971)  1527. 

8)  S.T.Picraux,  W.L. Brown  and  W.M. Gibson;  Phys. Rev.  B6  (1972)  1382. 

9)  K.Kinoshita,  T.Narusawa  and  W.M. Gibson;  Surface  Sci.  110  (1981)  33I. 

9a)  M.G.  Legally ,  in  Surface  Pr.vsics  of  Materials, voj. 2, ed. J.M.Bra&ely ,  Acaderdc  (1975  ! 

10)  D.P. Jackson,  3. M. Powell  and  G. Dolling;  Phys. Letters  51A  (1975)  67. 

11)  D.P. Jackson  and  J. n. Barrett ;Corput.  Phys.  Coznnun.  13.  (19T7 )  157. 

12)  P.M. Marcus,  J.E.Deruth  and  D.W.Jepsen;  Surface  Sci.  ^3  (1975)  501. 

12a)  R.E. Allen  and  F.W.de  Wette;  Phys.  Rev.  IBS  (i960)  1320. 

13)  J.E. Deruth,  P.M. Marcus  and  D.W.Jepsen;  Pays. Rev.  B15.  (1975)  1^60. 

13a)  P.E. Horten ,  R.L. Tapping  and  J. W.Goodale ;  Surface  Sci.  65  01977 )  13* 

lU)  P.H. Holloway  and  J.B. Hudson;  Surface  Sci.  0}  (197*0  123,  lUl. 

15)  R.G.Sreenk,  K.M.Tronp,  J.F.Van  der  Veen  and  F.W. Saris;  Surface  Sci. 95 


ptl9. 


(1980)  156. 


-29- 


Figure  Captions. 

Fig. I  Marble  model  of  bulk-like  Hi  (111)  surface  atomic  structure. 

Circles  of  different  sizes  correspond  to  atoms  in  different  planes. 

The  large,  medium,  and  snail  circles  in  the  top  viev  represent  the  atoms 
in  the  first,  second,  and  third  layers,  respectively.  The  side  viev  shows 
(llO)  planes  perpendicular  to  the  surface.  Note  the  incident  beam 
in  the  <!!!>  channeling  direction  sees  top  three  monolayers  whereas 
only  the  first  monolayer  in  the  <110:>  or  <10D>  direction.  A  and  3  in  the 
top  view  show  two  different  kinds  of  3-fold  symmetry  sites  for  oxygen 
adsorption. 

Fig. 2  Typical  aligned  backscattering  spectra  from  clean  and  oxygen  adsorbed 
Ni(ill)  surfaces. 

He+  ion  beam;  0.6  MeV ,  Q.h  pC  and  <110>  incidence.  The  exposure  to 
gas  is  100  L  (1  L  =  lxiO  ^  Torr*sec)  at  room  temperature  for  the  upper 
spectrum. 

Fig. 3  Angular  dependence  of  backscattering  spectrum  (a),  and  bulk  and  surface 


yields  (b). 

Tne  incident  He  ion  beam  is  at  1  KeV  and  in  the  vicinity  of  tr.e  <110> 
direction  in  the  {llO}  plane.  Tne  bulk  and  surface  yields  are  defined  as 
the  integrated  counts  in  the  segments  labelled  "B"  and  "S"  in  (a),  and  plot 
in  closed  and  open  circles  in  (b),  respectively.  Circled  numbers  in 
and  (b)  correspond  to  each  other  and  show  angular  netting  of  the  crystal 
with  resreet  to  the  incident  beam.  Data  were  taken  for  the  (►'r3x/3)R30‘' 


structure. 
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Fig. I  Energy  dependence  of  surface  peak  intensity  for  clean  Ki( 111 )-( lxl ) 
structure . 

(a):  <111>  incidence,  (b):  <110>  incidence.  The  uncertainty  of  absolute 
measurements  (+5  %)  are  shovn  by  error  bars.  Solid  curves  in  both  figures 
are  the  results  of  computer  simulations.  Simulation  (a)  assumes  a  bulk-like 
structure  and  isotropic  bulk-like  thermal  vibration,  i.e.  one-dimensional 
rms  amplitude  of  0.06S  X.  Simulation  (b)  assumes  a  bulk-like  structure 
and  an  enhanced  surface  thermal  vibration,  i.e.  one-dimensional  rms  amplitude 
of  0.081+  £  for  the  first  layer,  0.077  £  for  the  second  layer,  0.072  «  for 
the  third  layer,  and  0.C63  a  for  the  fourth  and  underlying  layers.  \ 

Fig. 5  Angular  dependence  of  <110>  surface  peak  intensity  at  1  KeV  for  clean 
Ni (lll)-(lxl )  structure. 

The  uncertainty  of  relative  measurements,  shovn  by  the  error  bars,  vas  estimated  to 
■be  ±3  £  and  ±6  £  near  the  centroid  and  the  peripheries  of  the  distribution,  respect i 

Solid  and  dashed  curves'  show  the  results  of  computer  simulations  in  which 

n 

a  bulk-like  structure  and  a  uniform  expansion  of  the  first  layer  by  0.05  « 
are  assumed,  respectively. 

Fig. 6  Variation  of  <111>  and  <110>  surface  peak  intensities  at  2  KeV  as  a  function 
of  exposure  to  oxygen.  Room  temperature  measurements. 

Exposures  which  give  most  intense  superstructure  diffraction  patterns  are 
inducated  by  arrows.  Error  bars  are  for  ±3  £  for  relative  SP  intensity 
measurements  and  ±10  £  for  oxygen  exposure  measurements . 

Fig. 7  Ehergy  dependence  of  <110>  surface  peak  intensity  for  Ni (lll)-(/3x/3)H3C?-0 
structure. 

Simulation  curves  assume  uniform  relaxation  of  the  first  monolayer  in  the 
direction  perpendicular  to  the  surface  by  various  amounts  indicated  in  the 
figure.  The  same  enhanced  thermal  vibration  model  as  in  Fig.U  is  assumed. 


>  \ 
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Fig. 8  Angular  dependence  of  <110>  surface  peak  intensity  at  1  KeV  for 
^  (a)  Ni(lll)-p(?x2)-0  structure  and  (b)  III ( 111 )-(  /Jx/F)R30°-0  structure. 

Simulation  curves  assume  the  enhanced  surface  thermal  vibration  (fig.h) 
and  uniform  expansion  of  the  first  monolayer  in  the  direction  perpendicular 
to  the  surface  by  various  amounts  indicated  in  the  figure.  The  exact  <110> 

\ 

direction  vas  determined  from  simultaneous  measurements  of  the  bulk  yield, 
a  Fig. 9  Angular  dependence  of  <110>  surface  peak  intensity  at  1  KeV  for  IJi(lll)- 

.  p(2x2)-0  structure  compared  to  simulations  of  non-uniform  surface  expansion 

models.  Simulation  (a)  assumes  the  expansion  of  all  nickel  atoms  in  the 
first  monolayer  by  0.15  *  (same  as  Fig. 8(a)).  Simulation  (b),  (c),  and  (fi) 
assume  that  3 A ,  1/2,  ar.d  l/b  of  the  nickel  atoms  in  the  first  monolayer 

D 

are  expanded  by  0.15  A,  respectively,  and  the  remainder  are  not  relaxed 
at  all.  The  enhanced  surface  thermal  vibration  is  assumed. 

Fig. 10  Atomic  structure  models  for  Ki(lll)-(lxl) »  -p(2x2)-0,  and  -( /Jx/3)R30°-0 
structures . 

Large  shadowed  circles  represent  oxygen  atoms  adsorbed  on  3-fold  symmetry 
sites  (site-B  in  Fig.l).  Small  closed  and  cpen  circles  represent  the  expanded 
and  ncn-expanded  nickel  atoms  in  the  first  monolayer,  respectively. 

The  second  and  underlying  layers  are  bulk-like  {not  shown  in  the  figure). 

Fig. 11  Energy  dependence  of  (a)  <111>  and  (b)  <11C>  surface  peak  intensities  for 
oxygen-saturated  Ki(lll)  surface  at  room  temperature.  Experimental  results 
for  the  clean  surface  are  also  shown  for  comparison.  Solid  curves  are 
just  visual  fits.  Exposure  to  oxygen :100  L. 


Fic-12  Variation  of  <110>  surface  peak  intensity  at  2  MeV  as  a  function  of 
exposure  to  oxygen.  420  K  (open  circles)  and  520  K  (closed  circles) 
measurements .  Increase  of  the'  SP  intensity  by  1  atom/rov  correspond 
to  displacements  of  one  monolayer  in  the  first-order  approximation 


for  <110>  incidence. 
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